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1 Introduction 



One of the main goals of the BABAR experiment is to study CP Violation in neutral B- 
mesons. The neutral B meson system is similar to the neutral Kaon system in that we have 
two flavor eigenstates that mix. However, the phenomenology is rather different. In the 
Kaon system the physics is driven by a large difference in the decay widths, such that there 
are two states, the and the Kg, that have substantially different lifetimes. In the Bd 
system the widths of the two states are very similar and instead the physics is dominated 
by the mass difference, Am^^, which controls the oscillation frequency of the B^B^ system. 

In the Standard Model, CP Violation arises from complex phases in the CKM matrix. By 
studying the interference between decays of B mesons that decays directly to a common final 
state and those that mix before they decay, we can study the phases of the CKM triangle[^. 

Unitarity of the CKM matrix allows the construction of unitarity triangles. Applying 
the unitarity constraint between the first and third generation gives the least degenerate 
triangle. The angles of this unitarity triangle are probed in CP Violation in B decays. 
Figure [1| shows the normalized CKM triangle, the studies of CP Violation in B —>■ J/tpKg 
allow us to determine sin 2/5. 

The decays of primary interest in this presentation are 5° —>■ J / ipKg and the very similar 
B^ —y ilj{2S)Kg. These decays have been dubbed the golden modes for measuring sin 2/5. 
These modes are experimentally fairly straight forward to construct, and have branching 
fractions that allow us to collect samples of events that are sufficiently large to study the 
time dependent asymmetry. But foremost, they are the golden modes due to the very small 
theoretical uncertainties as there are no penguin contributions with different weak phases. 

The time dependent rate for B J/ipK^ is given by 

/±(At;r, Ams,,sin2/5) = ^re"^!^*! [1 ± sin 2/3 x sinAms^At], (1) 

where +(— ) indicates that the -Btag, the other B from the T(4S') decay, was tagged as a 
B^ (B^). To experimentally fit the time distribution we need to account for two effects; 
finite resolution in the At determination and the possibility that the wrong tag was assigned 
to -Btag- The time resolution is handled by convoluting Eq. |I| with a resolution function, 
7^(At;d). The fraction of mistags, w, dilutes the sin 2/? measurement by P = 1 — 2w. 
Experimentally we will perform a fit to 

J'±(At;r, AmB^,Dsin2/3,a) = /±(At; P, Am^^, 1) sin 2/3) ® 7^(A^; a). (2) 

The time dependent CP asymmetry for the B^ J/ipK^ and B^ —>■ ilj{2S)Kg decays is 
given by 

where N^o (N^o) is the number of events where Bta.g was assigned a B^ (5°) tag. In these 
expressions At = tcp — ^tag refers to the time difference between the decay of Bcp and Btag- 
At the T(45') where the two B mesons are produced in a coherent state, the flavor of the B 
decaying to the CP eigenstate is determined by studying the flavor of the other B, as the two 



2 



B mesons are in a coherent P-wave state we know that the flavor of Bcp is of the opposite 
flavor of the B^ag at the time of the i?tag decay. Hence, At is the time Bcp evolved from 
when we knew its flavor. Further, at PEP-II the collisions are asymmetric, this allows us to 
measure the difference between the decay times of the two B mesons by simply measuring 
the separation in z of the two B decays. At ^ Az/c{P'y). 
The major steps in performing this analysis are 

• Reconstruct the CP eigenstates. 

• Measure the vertex resolution, 7l{At; a). 

• Determine vertex separation, which gives At 

• Tag the flavor of the other B, -Stag- 

• Measure the wrong tag fraction, w. 

• Perform likelihood fit to the At distribution to determine sin 2/3. 

As far as possible we try to determine the resolution function parameters and wrong tag 
fractions from data. 

2 The BABAR experiment 

The BABAR experiment is located at the PEP-II storage ring at SLAG. PEP-II collides 
electrons and positrons with energies of about 9.0 GeV/c^ and 3.1 GeV/c^ respectively at 
the center of mass energy of the T(4S') resonance. The produced T{4S) mesons have a boost 
of about j3'~f = 0.56. Since the first recorded collisions with the BABAR experiment on May 
26, 1999, PEP-II has produced excellent luminosity that have allowed BABAR to collect the 
worlds largest sample of B mesons at the T{4S). The analysis presented here is based on 9.0 
fb~^ recorded on the resonance and 0.8 fb~^ taken about 40 MeV/c^ below the resonance. 
This corresponds to about 10.5 x 10^ produced BB pairs. (The 2000 run of BABAR ended 
on October 30, 2000, and the total recoded luminosity in 2000 was 23 fb^^.) 

The BABAR experiment is described in detail elsewhere |@]. Here just a few key points 
of particular relevance to the measurement presented will be discussed. 

One of the key features of this experiment is that the produced B mesons have a boost 
of about /?7 = 0.56, and that the time difference between the B decays is measured by the 
separation in z position of the B decay vertices. The typical vertex separation between two B 
meson decays is 250 fim. The BABAR experiment has a 5 layered, double sided, silicon micro 
strip vertex detector capable of stand-alone tracking for low momentum particles, pt < 120 
MeV/c, not detected in the drift chamber. Figure |^ shows the measured track impact 
parameter resolution in ^ as a function of momentum. At higher momenta, where multiple 
scattering is negligible, the resolution is about 40 /im. For an exclusively reconstructed B 
the vertex resolution in z is typically about 40-60 /im. 
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Another unique feature of the BABAR experiment is the charged hadron identification 
system, the DIRC, Detection of Internally Reflected Cherenkov light. The DIRC detects 
Cherenkov photons that are produced in quartz bars and reflected out to a water tank 
instrumented with photo multipliers, see Figure |^. The DIRC has proved to work very 
well, we have achieved better than 3a K — n separation at momenta of 3 GeV/c. Further 
improvements are possible with a better understanding of the alignment of the DIRC with 
respect to the tracking system. 

Overall the BABAR experiment has performed very well, important for our ability to 
record data is our efficiency, typically BABAR is live 97% of the time when PEP delivers 
luminosity. Losses are due to background spikes that cause trips, ramp-up of voltages, and 
sporadic outages, e.g., of computing resources. After the data is recorded, it is typically 
reconstructed within 24 hours. 

2.1 Reconstruction of event samples 

The CP sample consists of events reconstructed in the following modes 

where the J/ip and the 'ip{2S) are reconstructed in both the e~^e~ and fi^^^ channels, the 
ipi^S) is also reconstructed in the J/iI)tt^tt~ channel. In the e^e~ decays, bremsstrahlung 
photon recovery is attempted. For more details about the event reconstruction see Ref. [^. 
We obtain a total of 168 events in these modes, as shown in Figures The yields and 
purities are summarized in Table |I]. 

To measure the performance of the tagging, and in particular to determine the wrong 
tag fraction, a sample of exclusively reconstructed B mesons is used. The modes used and 
the yields are summarized in Table ^. The reconstruction of these samples are detailed in 
Ref. § and §. 

2.2 Vertex separation 

For events where one exclusively reconstructed B meson has been found, the vertex of the 
other B decay is determined by trying to combine all other tracks in the event. Candidate 
tracks that form a good separate vertex, e.g., a K^, are combined to form a neutral candidate, 
which is used instead of the daughter tracks in the vertex determination. Tracks are removed 
if they contribute more than 6 to the of the vertex fit. Events are also rejected if \ Az\ > 3 
mm or if (Jaz > 400 fim. 

The time resolution function is parametrized as a sum of two Gaussians, 

7^(At;a) = ^-^exp(-(At-5,)V2a,^)- (4) 
i=i o-iV^vr 

The resolution parameters, (Tj, are taken as a scale factor. Si, times the calculated resolution 
based on the the tracking errors. The parameters for the second, wider, Gaussian is fixed 
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from Monte Carlo, the parameters for the first Gaussian is determined from the combined fit 
for mixing and the wrong tag fractions, see Section |2.3| . We also allow for a wide term, 
with a resolution of 1.8 ps, and no bias, to handle outliers. The parameters of the resolution 
function are summarized in Table ^. 



2.3 Flavor tagging 

Each event with an exclusively reconstructed decay is assigned a tag as a or a if 
the rest of the event satisfies the criteria for one of several tagging categories. These tagging 
categories are constructed such that each event will only belong to one category. The first 
category uses primary leptons to determine the fiavor. If the event contains an identified 
lepton, electron or muon, with center of mass momentum greater than 1.1 GeV/c the event 
is tagged as a B^ (5°) if the charge of the lepton is positive (negative). The second category 
uses charged particles identified as Kaons. If the sum of Kaons charges is positive (negative) 
the event is assigned a B^ (5°) tag. If the lepton tags and Kaon tags disagree no tag is 
assigned in these categories. The last two categories, NTl and NT2, are assigned based on 
the output of a neural network. The neural network combines information about Kaons, 
leptons, soft pions and the stiffest track in the event to form an output that distinguishes 
between B^ and B^ tags. The output from the neural network is shown in Figure |^ The 
two tagging categories are defined such that NTl corresponds to the events which has the 
best separation and NT2 to the events that has slightly worse separation. The events, in the 
middle of Figure ^ which have very little separation are not used. 

The figure of merit for each tagging category is the effective tagging efficiency, Qi = 
ej(l — 2wiY = ^i'^h where is the fraction of events assigned to category i and Wi is the 
fraction that had the wrong tag assigned. 

To determine the wrong tag fraction we use the sample of exclusively reconstructed 
hadronic B meson decays. These decays tag the flavor of the decaying B, so by performing 
a combined tagging and mixing analysis we can determine the wrong tag fraction for each 
category. For the events in each tagging category we perform a flt to 

n±{At;r,AmB^,V,a) = h:e~^^^^^[l ±V x cosAm^^At] ® 7^(A^;a), (5) 
where + are unmixed events and — are mixed events. A log-likelihood is formed by 

\nC = Y,[J2 ^nn+{At;T,AmB„V„a)+ \nn-{At;T,AmB„V„a)] (6) 

i unmixed mixed 

where i runs over the tagging categories. Additional terms are added to the probability 
density functions to describe the contributions from backgrounds, details are given in . 

The results of the flt are shown in Table ^. The total tagging efficiency is 76.7 ± 0.5% 
with an effective tagging efficiency, Q, of 27.9 ± 0.5% (statistical errors only). 

When the tagging algorithm is applied to the sample of 168 CP events 120 events were 
assigned a flavor tag. Table ^ shows a break down per mode and per tagging category of the 
events in the CP sample. Of the 120 events tagged, 70 were B^ and 50 were B^ tags. 
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2.4 Systematics 

Systematic errors were considered from many different sources; input parameters to the 
likelihood fit, uncertainties in the time resolution function and wrong tag fractions. The 
lifetime was fixed to the PDG [^] value r^o = 1.548 ps, AniB^ was also fixed to the PDG 
value, ^tTiB^ = 0.472 hps^^. Varying the values of these parameters give the uncertainties 
on sin 2/3 for r^o and Am^^ of 0.002 and 0.015 respectively. 

The time resolution function was determined in a high statistics sample of fully recon- 
structed events. We vary the parameters in the time resolution function by 1 statistical 
standard deviation and assign a systematic error on sin 2/5 of 0.019. To study the sensi- 
tivity on the bias in At, we allowed the bias of the second Gaussian to increase to 0.5 ps. 
This results in a change of 0.047 on sin 2/3 and is assigned as a systematic uncertainty. The 
sensitivity to this bias is due to the different number of events that are tagged as B'^ and 

The mistag fractions are determined from exclusively reconstructed 5° and 5° mesons. 
Sources of systematic uncertainty come from the presence of backgrounds in these samples 
and possible differences between the tagging performance in the CP sample and the hadronic 
samples. The details about the accounting of backgrounds in the hadronic samples are given 
in Ref. [^]. The systematic uncertainty on sin 2/? from the measured mistag fractions is 
estimated to be 0.053. A rather conservative systematic error of 0.050 on sin 2/3 is assigned 
for a possible difference between the tagging performance between the CP sample and the 
exclusive hadronic sample. 

The CP sample is estimated to contain a background fraction of (5 ± 3)%. Backgrounds 
from, e.g., u, d, and s continuum events contributes primarily at small values of At and 
hence do not contribute much to the determination of sin 2/3. We estimate that the effective 
background is 3% and correct for the background by increasing the apparent asymmetry 
by 1.03. A fractional systematic error of 3% is assigned on the asymmetry to cover the 
uncertainty in the size of the background as well as any possible CP asymmetry in the 
background. 

The systematic errors are summarized in Table |^, and a total systematic uncertainty of 
0.09 on sin 2/3 is obtained. 

2.5 Results 

The analysis was carried out blind to eliminate any possible experimenters bias. The blinding 
technique hid both the result of the likelihood fit for sin 2/3, as well as the CP asymmetry in 
the At distribution. The error on the asymmetry was not hidden. The value of sin 2/3 was 
hidden by adding an arbitrary offset and flipping the sign. The CP asymmetry in the At 
distribution was hidden by adding an offset and, on an event by event basis, multiply At 
with the sign of the tag. 

This allowed us to carry out many systematic studies while keeping the value of sin 2/3 
hidden. In particular, the whole analysis procedure, including event selection, was fixed 
prior to unblinding the value of sin 2/3. 

Using Eq. we perform a likelihood fit to determine one single parameter, sin 2/3. The 
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mistag fractions, Wi, and the resolution function parameters, a are taken from Tables ^ and ^ 
respectively. We obtain 

sin 2/3 = 0.12 ±0.37 ±0.09. (7) 

Table |^ shows a breakdown of the fit for different modes and tagging categories. Figure ^ 
shows the At distributions for and tagged events, the different yields of B'^ and B^ 
tags is apparent in this plot as well as in Figure |^, which shows the raw asymmetry. 

To validate the analysis several cross checks have been made. In particular, we have 
used the charmonium non-CP samples such as B^ J/ipK^ and the self tagging B^ mode, 
J/ipK*^, K*^ — > K^TT^, to establish that we do not observe any time dependent asymmetry 
in these modes. We also use the large samples of exclusively reconstructed charged and 
neutral B mesons in this check. Table ^ summarizes the fits for an apparent CP asymmetry 
in these modes. The result is consistent with no time dependent asymmetry in these modes. 

Other important validation measurements are the lifetime and mixing measurements 
^, also discused in detail in Rainer Bartoldus contributions to this conference. 

Figure [TO shows the constraints in the p — rj plane and the preliminary BABAR mea- 
surement of sin 2/3. The allowed area in the plot does not include any constraint from our 
measurement of sin 2/3. 



3 Conclusion 

BABAR has reported preliminary results based on a dataset of 9.0 fb^^ recorded at the 
T{4S) resonance from January to July in 2000. The analysis used the modes B^ J/ipKg 
and B^ ip{2S)K^ and studied the time dependent asymmetry in events where the flavor 
of the other B meson was determined. The preliminary result 

sin 2/3 = 0.12 ±0.37 ±0.09 (8) 

was obtained. An analysis of the full 2000 dataset is now in progress, other modes, including 
B^ — > J/ipKi^, is to be included to further improve the precision of the measurement of 
sin 2/3. Data collected over the next few years should significantly improve the precision of 
this measurement. BABAR and PEP-II will start running again in February of 2001. 
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Final state 


Yield 


Purity (%) 


J/i)Ks {Ks TT+TT-) 


121 


96 




19 


91 


^{2S)Ks {Ks TT+TT") 


28 


93 



Table 1: The yield of events in the different modes for the CP sample. 



Final State 


Yield 


Purity (%) 


D*-7r+ 


622 ± 27 


90 


D*-p+ 


419 ±25 


84 


D*-at 


239 ± 19 


79 


D-TI+ 


630 ± 26 


90 




315 ± 20 


84 


D-ai 


225 ± 20 


71 


Total 


2438 ± 57 


85 




1755 ±47 


88 


D*07T+ 


543 ± 27 


89 


Total 


2293 ± 54 


88 


D*-i+iy 


7517 ± 104 


84 



Table 2: The yield of events in the different modes for B mesons reconstructed in hadronic 
and semileptonic modes. 



9 



Tagging Category 


ei%) 


w (%) 


Qi%) 


Lepton 


11.2 ±0.5 


9.6 ± 1.7 ± 1.3 


7.3 ± 0.3 


Kaon 


36.7 ±0.9 


19.7 ± 1.3 ± 1.1 


13.5 ±0.3 


NTl 


11.7±0.5 


16.7 ±2.2 ±2.0 


5.2 ±0.2 


NT2 


16.6 ±0.6 


33.1 ±2.1 ±2.1 


1.9 ±0.1 


all 


76.7 ±0.5 




27.9 ±0.5 



Table 3: The efficiency and wrong tag fractions in the different tagging categories as 
determined from the combined tagging and mixing fit to the hadronic and semileptonic 
event sample. 









CP sample 




— s> 71 


+7r-) 




— s> n 






— > n 


+71-) 


(tagge 


d) 






all 




B^ 


all 


B^ 


B' 


all 






all 


Electron 


1 


3 


4 


1 





1 


1 


2 


3 


3 


5 


8 


Muon 


1 


3 


4 











2 





2 


3 


3 


6 


Kaon 


29 


18 


47 


2 


2 


4 


5 


7 


12 


36 


27 


63 


NTl 


9 


2 


11 


1 





1 


2 





2 


12 


2 


14 


NT2 


10 


9 


19 


3 


3 


6 


3 


1 


4 


16 


13 


29 


Total 


50 


35 


85 


7 


5 


12 


13 


10 


23 


70 


50 


120 



Table 4: The result of tagging applied to the CP sample and broken down by tagging 
category and mode. 



Parameter 


Value 


6i (ps) 


-0.20 ±0.06 


from fit 




1.33 ±0.14 


from fit 


L (%) 


1.6 ±0.6 


from fit 


/2 (%) 


25 


fixed 


S2 (ps) 





fixed 


^2 


2.1 


fixed 



Table 5: The parameters of the resolution function as determined from the sample of fully 
reconstructed hadronic B° candidates. 
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bclllipKJ 


Sill zp 


Ur sample 


0.12±0.37 


J /vjKs [-Us TT TT j events 


— U.iU ± [jAZ 


otlior CP cn'Oiits 


0.87 ± 0.81 


Lepton 


1.6 ± 1.0 


Kaon 


0.14 ±0.47 


NTl 


-0.59 ±0.87 


NT2 


-0.96 ±1.30 



Table 6: The result of the fit for sin 2/? broken down by event and tagging category. 



Systematic Error 


Uncertainty on sin 2/3 




0.002 




0.015 


Az resolution 


0.019 


time resolution bias 


0.047 


measured mistag fraction 


0.053 


CP versus non-CP sample 




mistag fraction 


0.050 


versus mistag fraction 


0.005 


background in CP sample 


0.015 


total systematic error 


0.091 



Table 7: Summary of sources of systematic errors in the determination of sin 2/3. 



Sample 


Api)areiit CP as}'iiiiiietry 


Hadronic charged B decays 


0.03 ±0.07 


Hadronic neutral B decays 


-0.01 ±0.08 




0.13 ±0.14 




0.49 ± 0.26 



Table 8: Results of fitting for apparent CP asymmetries in various charged or neutral flavor 
tagging B samples. 
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Figure 1: The CKM triangle formed by unitarity between the first and third generation. 
The study of time dependent asymmetry in J / tjjKg allows us to determine sin 2/3. 
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Figure 2: The impact parameter resolution in z for the silicon vertex tracker as measured 
by cosmic muons. 
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1—91 mm^ MOmm 
"4 5m 1 j. 1.17m 



4 X 1 .225 m Bars 
glued end-to-end 

Figure 3: The principle of the DIRC; incident charged particles emits Cherenkov radiation 
as they pass through the quarts bars. The Cherenkov photons are reflected to the end of the 
quartz bar and out in a water tank, which is instrumented with photomultipliers to detect 
the light. 
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Figure 5: The signal in the i?^ 
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Figure 6: The signal in the —>■ ip{2S)Kg, Kg tt+tt mode. 
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P(B") 



Figure 7: The output from the neural network for B flavor tagging. The data is shown as 
points with error bars and the histogram shows the Monte Carlo. Events close to 1 (0) are 
likely to be (5") tags. Events near 0.5 provide no taging information and are not used. 
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Figure 8: The At distributions for events in the CP sample tagged as -B" and B^ respectively. 
The asymmetry in the number of 5° and S° tags are apparent in the plot. The shift to 
negative values of At is due to the small bias, about 0.2 ps, in the time resolution. 
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Figure 9: The raw — asymmetry ( A^^o — A^^o ) / ( A^^o + A^^o ) , with binomial errors, as 
a function of At. The sohd hue is the asymmetry for our central value of sin 2/3. The dotted 
lines represents one statistical standard deviations from the central value. 
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Figure 10: The constraints in the p — fj planes, and the result of our determination of 
sin 2/? overlayed. The following measurements are used: \Vcb\ = 0.0402 ± 0.017, |V^b/V^fe| = 
{\Vub/Vcb\) ± 0.0079, Ams, = 0.472 ± 0.017 ^s'^ and |ex| = (2.271 ± 0.017) x lO'^, and 
for Am^^ the set of amplitudes corresponding to a 95%CL limit of 14.6 /ips^^. We scan the 
model-dependent parameters {\Vuh/Vch\), Bk, fs^ \l Bb^ and in the range [0.070, 0.100], 
[0.720, 0.980], [185, 255] MeV and [1.07, 1.21], respectively. sin2/? = 0.12 ± 0.37(stat) 
is represented by cross-hatched regions corresponding to one and two statistical standard 
deviations. 
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